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ABSTRACT 


The  reduction  and  conversion  of  airborne  gravity  gradi ometer  data  to  ground 
level  estimates  of  the  gravity  disturbance  vector  is  currently  of  considerable 
interest  in  support  of  short  wavelength  gravity  modeling.  A  pressing  problem 
is  the  need  for  an  accurate  procedure  for  the  downward  continuation  of  data 
acquired  at  altitude  by  the  Airborne  Gravity  Gradiometer  Survey  System  (GGSS). 
As  part  of  ongoing  investigations,  a  prism  method  has  been  used  to  calculate 
the  effect  of  topography  on  the  gravity  disturbance  vector  and  the  five  inde¬ 
pendent  second-order  gravity  gradients.  Calculations  of  the  contribution  of 
topography  to  the  magnitudes  of  these  gravimetric  parameters  were  made  at  both 
surface  and  elevated  points  in  the  Wichita  Mountains  of  Oklahoma.  Computations 
were  made  utilizing  Digital  Terrain  Elevation  Data  (DTEu )  with  an  assumed  con¬ 
stant  density  of  2.67  grams/centimeters'3  for  the  topographic  masses.  Results 
are  presented  which  reflect  the  use  of  DTED  sets  of  different  horizontal  extent 
and  grid  interval. 
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THE  EFFECT  OF  TUPOGRAPHY  ON  AlkBORNE  GRAVITY  l, RADIOMETER  DATA 


I.  INTRODUCTION 

The  field  testing  of  the  airborne  Gravity  Gradioineter  Survey  System  (GDSSl, 
being  built  by  Bell  Aerospace/Textron  for  the  Defense  Mapping  Agency,  is 
scheduled  to  commence  before  midyear.  Therefore,  methods  for  validation  or 
the  system's  ability  to  map  the  local  gravity  field  to  high  detail  through 
gravity  gradient  measurements  and  their  subsequent  downward  continuation/ 
conversion  into  surface  gravity  disturbance  components  are  or  pressing 
interest.  Validation  can  in-part  be  accomplished  by  estimating  the 
influence  of  local  topography  on  the  radial  disturbance,  the  deflection 
components,  and  the  second-order  gravity  gradients  at  both  surface  and 
aloft  stations.  This  report  summarizes  the  results  of  computing  topo¬ 
graphic  terrain  effects  from  the  topography  above  mean  sea  level.  The 
terrain  effects  are  calculated  on  the  basis  of  homogeneous  rectangular 
prisms  which  model  the  terrain  masses  with  an  assumed  constant  density. 

The  surface  computation  points  coincide  with  the  two  astro-geodetic  sta¬ 
tions  located  in  the  Wichita  Mountains  of  Oklahoma.  One  is  near  Sunset 
Peak  and  the  other  is  on  Mount  Scott  as  shown  in  Figure  1.  Terrain  effects 
were  computed  at  thesp  surface  stations  and  at  points  directly  overhead  at 
altitude  5000  feet  above  the  geoid.  Digital  Terrain  Elevation  Data  (DTEDI 
supplied  the  topographic  model  needed  to  compute  the  terrain  effects.  The 
main  goal  of  the  study  was  to  determine  which  DTEu  field  should  make  tin  an 
inner  grid  zone  and  to  what  radial  extent,  outward.  There  was  also  a  need 
to  establish  the  coarsest  DTEU  representa i on  permissible  for  the  outer  zone 
and  its  span  of  coverage  for  adequate  modeling  of  the  local  terrain  effects 
on  selective  gravimetric  quantities. 

II.  DISCUSSION 
a.  Objective 

A  major  part  of  the  short  wavelength  variation  of  the  qravity  field  in  a 
local  area  is  due  to  topography.  Therefore,  especially  in  mountainous 
areas,  one  would  anticipate  the  nepd  for  the  bpst  available  topographic 
elevation  data  in  the  immediate  vicinity  of  a  computation  point.  At  some 
further  distance  beyond  this,  a  coarser  terrain  representation  could  he 
utilized  to  make  the  computational  process  more  efficient  with  a  minimal 
affect  on  results.  With  this  premise,  terrain  effects  were  computed  for 
the  stations  shown  in  Figure  1  to  establish  the  radius  of  an  inner  zone 
for  3"  point  DTEU  which  is  our  finest  grain  DTEU.  The  computations  also 
allow  the  determination  of  the  largest  mean  UTED  representation  that  is 
permissible  for  an  outer  zone  and  its  span  of  coverage  for  adequate 
modeling  of  the  local  terrain  effects  on  selective  gravimetric  quantities. 
These  quantities  include  the  radial  gravity  disturbance,  the  deflection 
of  the  vertical  components,  and  the  second-order  gravity  gradients. 

Inner  zone  modeling  with  mean  DTEU  will  indicate  whether  there  is  a  need 
for  the  exclusive  use  of  point  DTEU  for  this  area.  Insiqht  into  upward 
continuation  effects  is  afforded  with  the  inclusion  of  aloft  computation 
points. 
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b.  Method  of  Analysis 

A  modified  version  of  the  Kene  Forsherg  prism  program  was  run  to  compute 
the  topographic  terrain  effects.  The  program  is  discussed  in  keferenc.*  2. 
The  right-handed  coordinate  system  was  centered  at  each  computation  sta¬ 
tion  with  X  pointing  east,  Y  pointing  north,  and  Z  pointing  up.  Prisms 
were  formed  from  a  geoid  base  up  to  the  DTEU  topography  of  an  assumed 
density  of  2.67  gm/cnr  as  illustrated  in  Figure  2.  The  integration  was 
performed  numerically  using  the  distances  from  station  coordinates  to  the 
eight  corners  of  each  prism.  The  surface  station  height  was  offset  1  cm 
to  avoid  the  central  prism  from  being  automatically  eliminated  from  the 
program  and  to  avoid  the  station  from  being  located  within  the  central 
prism  boundaries. 

c.  UTEU  Fields  Used 

A  3"  point  DTEU  field  was  built  around  the  stations  from  the  same  data 
that  fJMAAC  sent  to  Bell  Aerospace/Textron.  Also  created  were  4”,  12". 

1 5 “ ,  30",  1',  and  3‘  mean  fields  from  the  3“  point  data.  When  using 
both  an  inner  field  and  an  outer  field  around  a  station,  the  prisms  must 
properly  fit  together.  If  the  inner  field  is  the  3'1  UTEU,  then  the  outer 
mean  field  must  be  a  multiple  of  6"  for  the  prisms  to  fit  together. 

d.  Formulation 

Terrain  mass  modeling  is  accomplished  with  homogeneous  rectangular  prisms 
in  the  calculations.  Gravitational  formulas  for  such  prisms  are  known 
from  MacMillan's  work  on  potential  theory  in  Reference  1.  The  prism 
dimensions  for  this  study  were  controlled  by  the  description  of  the  topo¬ 
graphy.  Figure  3  shows  the  indefinite  integral  solutions  used  to  compute 
the  exact  evaluation  of  ten  gravimetric  guantities  at.  the  computation 
point  P  in  Figure  2  for  a  single  prism.  When  the  separation  between  a 
prism  and  the  computation  point  permits,  approximate  prism  formulations 
are  used  instead  nf  the  exact  ones.  Figure  4  gives  details  on  how  a 
series  expansion  for  a  prism's  potential  may  be  derived  by  formulating 
the  reciprocal  distance  "r"  as  an  infinite  series  in  terms  of  the  Legendre 
polynomials.  Using  only  the  first  few  terms  of  the  potential  series,  all 
first  and  second-order  gravity  gradients  may  be  found  by  simple  differen¬ 
tiation.  The  exact  formulation  is  normally  used  tor  the  central  prism 
while  the  approximate  formulation  is  used  for  all  remaining  prisms  to 
obtain  the  desired  gradients. 


III.  RESULTS 

a.  Tabular  Output 

Tables  1  and  2  show  the  effects  of  topqraphic  modeling  by  different  UTFl) 
fields  within  a  12'  radius  about  the  two  surface  stations.  The  various 
terrain  fields  used  were  3",  9",  12",  15",  30",  1' ,  and  3'  UTEu.  The 
same  investigation  is  repeated  in  Tables  3  and  4  tor  the  computation 
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points  at  elevation  5000  feet  above  the  geoid  and  directly  overhead. 
Terrain  effect  variations  due  to  unit  step  increases  over  the  interval. 

6'  to  13',  in  the  radius  of  a  central  zone  of  3"  point  UTEU  are  shown  in 
Tables  5  and  6  for  the  ground  level  stations.  In  the  remaining  tables, 
terrain  effects  are  accounted  for  by  employing  the  3"  point  UTEU  in  an 
inner  zone  and  one  of  the  mean  DTEU  fields  in  an  outer  zone.  It  was 
discovered  that  the  prism  program  requires  the  outer  grid  field  to  he  a 
multiple  of  6"  for  the  prisms  to  properly  fit  together.  This  is  the  same 
as  saying  that  if  one  extends  the  31'  field,  then  the  outer  mean  field 
center  coordinates  must  coincide  with  a  3"  point.  Examples  of  outer 
fields  would  then  be  12“,  30",  1',  and  3‘. 

Tables  7  and  8  indicate  variation  in  terrain  effects  caused  by  mean  UTEU 
representations  in  an  outer  zone  of  radius  30'  and  an  inner  zone  radius 
of  12*.  The  means  considered  for  the  outer  zone  terrain  modeling  were 
12",  30",  1'  and  3'  LlTEQ.  Each  table  displays  results  for  one  of  the 
surface  points.  Sunset  Peak  or  Mount  Scott,  and  the  related  overhead 
point.  Tables  9  and  10  reflect  the  usage  of  1'  mean  UTEU  in  two  expanded 
outer  zones,  one  of  radius  2°  and  the  other  2.5°  with  an  inner  zone  radius 
of  12'.  These  are  again  composite  tables  for  a  pair  of  ground/aloft 
stations.  For  the  ground  level  stations.  Tables  11  and  12  indicate  ter¬ 
rain  variations  due  to  inner  zone  modeling  with  different  radial  extents 
over  the  interval  1*  to  10'  while  the  outer  1'  UTEU  zone  was  extended  out 
to  a  radius  of  2.5°.  Tables  13  and  14  are  composite  tables  for  a  pair 
of  ground/aloft  stations  where  the  terrain  modeling  was  accomplished  hy 
an  inner  region  of  3“  point  UTE0  and  an  outer  zone  of  3'  mean  UTEU.  The 
radii  of  the  inner  and  outer  regions  were  5'  and  2.5°,  respectively,  in 
one  case  while  6'  and  3°  in  another. 

b.  Analysis 

From  Tables  1  through  4,  it  is  obvious  that  the  Sunset  Peak  station  is 
not  on  a  very  sharp  peak  and  that  the  Mount  Scott  station  is  on  a  vptv 
sharp  peak  which  changed  by  200  meters  when  using  a  1'  mean  elevation 

field.  Tables  1  through  4  make  it  obvious  that  the  inner  grid  must  he 

the  3"  point  UTEU.  Tabies  5  and  6  show  that  the  maior  contribution  to  the 
terrain  effects  is  located  very  close  to  the  stations.  Tables  7  and  8 
show  that  l'  and  3'  mean  DTEU  outer  qrids  yield  the  same  results  as  a  1?" 

mean  outer  grid.  Tables  9  and  10  show  that  an  outer  radius  of  2.5°  allows 

the  final  effects  to  converge  within  acceptable  limits.  Tables  11  and  12 
show  that  the  inner  grid  radius  may  be  reduced  to  at  least  5'  and  still 
yield  very  acceptable  results  with  a  2.5°  outer  radius.  Tables  13  and  14 
show  that  an  outer  grid  of  3'  mean  UTEU  yields  acceptable  results  when 
compared  with  the  1'  mean  outer  grid  results. 

When  an  outer  grid  is  used,  the  inner  grid  values  are  extended  further 
than  they  were  in  Tables  1  through  4.  For  instance  when  the  3'  tield  was 
the  outer  grid  and  the  inner  radius  was  5'  as  in  Table  13,  the  portion  of 
fine  grid  used  in  computations  agrees  with  that  of  the  8'  radius  of  Table 
5.  Another  example  is  when  the  1'  field  was  the  outer  grid  and  the  inner 
radius  was  6'  as  in  Table  11.  This  agrees  with  the  portion  of  fine  grid 
used  in  computations  for  the  7'  radius  of  Table  5.  The  inner  field  is 
then  extended  by  one  unit  of  the  outer  grid  spacing  to  form  prisms  that 


fit  next  to  the  prisms  from  the  outer  field.  This  explains  why  the  outer 
grid  must  be  a  multiple  of  fi“  for  the  prisms  to  fit  together.  This  also 
explains  the  increase  in  number  of  prisms  for  Tables  7  and  8.  The  3" 
field  was  extended  further  out  by  1'  and  3'  which  caused  an  increase  in 
number  of  prisms.  However,  as  one  increases  the  outer  radius  to  3°,  a 
decrease  in  number  of  prisms  would  be  seen  by  usinq  a  3'  outer  mean  rield 
rather  than  an  outer  field  of  higher  density.  Computer  time  is  thus 
reduced  by  using  the  3‘  outer  field. 


IV.  CONCLUSIONS  AND  FUTURE  CONSIDERATIONS 

The  inner  field  must  be  the  3"  point  OTEU  with  an  inner  radius  of  at  least 
5‘.  The  outer  field  may  be  3‘  mean  DTEL)  with  an  outer  radius  of  at  least 
2.5°.  If  the  field  is  big  enough  it  would  be  advisable  to  use  an  outer 
radius  of  3°  and  an  inner  radius  of  6'  in  qeneral .  This  would  allow 
convergence  of  the  topographic  terrain  effects  within  acceptable  limits. 
The  3“  inner  grid  contributes  the  ma.ior  portions  ot  the  gradient  values 
which  means  that  the  3"  elevations  near  the  station  coordinates  must  he 
as  accurate  as  possible.  An  inner  grid  that  is  a  less  dense  mean  field 
would  not  yield  the  proper  elevations  near  the  station.  The  outer  qnd 
may  be  thought  of  as  fine-tuning  the  values  until  convergence  is  achieved 
within  acceptable  limits.  The  3"  DTED  tapes  are  written  in  3'  blocks 
which  makes  it  convenient  to  form  3'  mean  UTEU  fields  for  the  outer  grid. 

One  must  consider  the  computer  time  saved  by  using  3‘  mean  UTEU  fields 
for  the  outer  grid  as  large  numbers  of  stations  are  computed.  The  number 
of  prisms  is  cut  down  without  sacrificing  accuracy.  About  90  seconds  of 
CPU  time  is  added  per  station  with  3'  mean  UTEU  out  to  a  radius  of  3°  and 
3"  point  UTEU  out  to  a  radius  of  6'.  The  radial  gravity  disturbance  and 
second-order  gradients  converged  within  desired  limits.  The  deflection 
components  will  not  converge  as  the  outer  grid  radial  extent  is  increased. 
Conversion  of  the  deflections  into  the  other  two  disturbance  components 
does  not  qive  near  the  magnitude  of  the  radial  disturbance.  He  see  the 
high  frequency  nature  of  the  second-order  gradients  by  their  rapid  con¬ 
vergence.  The  much  lower  frequency  nature  of  the  deflection  components 
will  not  allow  convergence  in  the  computation  of  the  terrain  effects. 

In  the  future  it  would  be  desirable  to  augment  the  analysis  with  more 
stations  in  the  Wichita  mountains.  The  harmonic  content  ot  topographic 
terrain  effects  needs  to  he  determined.  Terrain  effect  computation  hy 
the  FFT  method  may  be  investigated  in  the  future.  Comparisons  of  actual 
gradiometry  data  against  simulated  gradients  derived  from  UTEU  fields  will 
need  to  be  made.  It  would  be  desirable  to  incorporate  error  propagation 
into  the  topographic  terrain  effects  programs  as  errors  in  the  UTEU  and 
density  assumption  become  better  defined.  A  maior  goal  is  to  he  able  to 
add  harmonic  quantities  to  terrain  effects  for  prediction  of  gravity  in 
inaccessahle  areas. 
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MU U NT  SCOTT  Et FECTS  USING  ONLY  3"  INNER  FIELD 
AND  STATION  ELEVATION  LOCATED  ON  TERRAIN  SURFACE 
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MOUNT  SCOTT  EFFECTS  USING  A  3"  INNER  FIELD  OF  12' 
RADIUS  AND  A  VARIABLE  OUTER  FIELD  WITH  A  RADIUS  OF 
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